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Abstract 

The multiple reference state structure of the Z„ Belavin model with non-diagonal 
boundary terms is discovered. It is found that there exist n reference states, each of 
them yields a set of eigenvalues and Bethe Ansatz equations of the transfer matrix. 
These n sets of eigenvalues together constitute the complete spectrum of the model. In 
the quasi-classic limit, they give the complete spectrum of the corresponding Gaudin 
model. 
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1 Introduction 



Our understanding of quantum phase transitions and critical phenomena has been greatly 
enhanced by the study of exactly solvable models (integrable models) pQ. Such exact results 
provide valuable insights into the key theoretical development of universality classes in areas 
ranging from modern condensed physics [2] to string and super-symmetric Yang-Mills the- 
ories [3]. Among solvable models, elliptic ones stand out as a particularly important class 
due to the fact that most other ones can be obtained from them by some trigonometric or 
rational limits. In this paper, we focus on the elliptic Z n Belavin model [I] with integrable 
boundary conditions, with the celebrated XYZ spin chain as the special case of n = 2. 

Two-dimensional integrable models have traditionally been solved by imposing periodic 
boundary conditions. For such bulk systems, the quantum Yang-Baxter equation (QYBE) 

R l2 {u x - u 2 )Rxz{u 1 - u 3 )R 23 {u 2 - u 3 ) = R 23 {u 2 - u 3 )R 13 (ui - u 3 )R 12 {ui - u 2 ), (1.1) 

leads to families of commuting row-to-row transfer matrix which may be diagonalized by the 
quantum inverse scattering method (QISM) (or algebraic Bethe Ansatz) [5]. 

Not all boundary conditions are compatible with integrability in the bulk. The bulk 
integrability is only preserved when one imposes certain boundary conditions. In [6J, Sklyanin 
developed the boundary QISM, which may be used to describe integrable systems on a finite 
interval with independent boundary conditions at each end. This boundary QISM uses a 
new algebraic structure, the reflection equation (RE) algebra. The solutions to the RE 
and its dual are called boundary K-matrices which in turn give rise to boundary conditions 
compatible with the integrability of the bulk model [6j [7J [8]. 

The boundary QISM has been successfully applied to diagonalize the double-row transfer 
matrices of various integrable models with non-trivial boundary conditions mostly corre- 
sponding to the diagonal K-matrices. The problem of diagonalizing the double-row transfer 
matrix for most general non-diagonal K-matrices by the algebraic Bethe Ansatz has been 
long-standing due to the difficulty of finding suitable reference states (or pseudo-vacuum 
states). Recently, much progress has been made for the open XXZ spin chain. Bethe Ansatz 
solutions for non-diagonal boundary terms where the boundary parameters obey some con- 
straints have been proposed by various approaches^ [9l [10l [TTJ [TJl [131 El [E3 OS] • It has been 

1 Solutions with arbitrary boundary parameters were recently proposed by functional Bethe Ansatz [17] 
and q-Onsager algebra [18) . However, it seems highly non-trivial to rederive these results in the framework 
of algebraic Bethe Ansatz . 
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found that in order to obtain the complete spectrum of the model two sets of Bethe Ansatz 
equations and consequently two sets of eigenvalues are needed [HI [15] , in contrast with the 
diagonal boundary case [6]. This suggests that in the framework of algebraic Bethe Ansatz 
there should exist two reference states corresponding to the two sets of Bethe Ansatz equa- 
tions and eigenvalues. Such multiple reference state structure was confirmed in our recent 
work [20]. 

However, for models related to higher rank algebras [211 1221 [231 [2H (23 [2S] only one 
reference state, and consequently only one set of eigenvalues and Bethe Ansatz equations 
of their transfer matrices, have been constructed so far. It is natural to expect that there 
exist extra reference states, giving rise to extra sets of eigenvalues and associated Bethe 
Ansatz equations for such models. In this paper we investigate the multiple reference state 
structure for an elliptic model related to algebra - the Z n Belavin model with general 
non-diagonal boundary terms. It is found that there actually exist n reference states for 
such a model. Each of these reference states yields a set of eigenvalues and corresponding 
Bethe Ansatz equations of the transfer matrix. The n sets of eigenvalues together constitute 
the complete spectrum of the model. In the quasi-classic limit, they give the corresponding 
spectrum of the associated Gaudin model [27] . 

The paper is organised as follows. In section 2, we briefly review the Z n Belavin model 
with integrable open boundary conditions, which also serves as an introduction to our notion 
and basic ingredients. In section 3, we introduce the intertwiner vectors and the correspond- 
ing face-vertex correspondence relations which will play key roles in transforming the model 
in "vertex picture" to the one in the "face picture" . In section 4, after finding the n reference 
states, we use the algebraic Bethe Ansatz to obtain the corresponding n sets of eigenvalues 
and the associated n sets of Bethe Ansatz equations of the transfer matrix of the model. In 
section 5, we take the quasi-classic limit to extract the spectrum of the associated Gaudin 
operators. Section 6 is for conclusion. The Appendix provides the definitions of some ele- 
mentary functions appeared in section 4 and 5. 
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2 7j n Belavin model with open boundaries 

Let us fix a positive integer n > 2, a complex number r such that Im{r) > and a generic 
complex number w. Introduce the following elliptic functions 
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C(«) = ^{\na(u)}. 

Among them the a-functiorj^ satisfies the following identity: 

a{u + x)a{u — x)a{v + y)a{v — y) — cr(u + y)a{u — y)&(v + x)cr(v — x) 
= cr(u + v)a{u — i>)cr(:r + y)a(x — y). 

Let g, h, be n x n matrices with the elements 

hij = S i+ ij, gij = u'dij, with u = e « , Z n . 

For any a = (ai, 02), oc\, 0,2 G Z n , one introduces annxn matrix J a by 



g a2 h a \ 



and an elliptic function cr a (ii) by 

cr a (ii) = 6* 
The Z n Belavin R-matrix is [U [28] 



1 I m 1 

2 ~ n 

1 I Q2 



(w,r), and (T(o,o)(u) = a(u) 



R B (u) 



a(w) 



a a (u+ J2) , 



cr(w + w) ^— < ncTQ,( H ) 



The R-matrix satisfies the QYBE (II. ip and the properties [28 
Unitarity : 

Crossing-unitarity : (R B )\\(-u - nw)(R B )%(u) 



R B 2 (u)Rg(-u) =id, 



<j(u)a(u + nw 



a(u + w)a(u + niy — w) 
Quasi-classical property : R B 2 (u)\w^o = id. 



(2.1) 

(2.2) 
(2.3) 



(2.4) 



(2.5) 
id, (2.6) 
(2.7) 



2 Our cr-function is the ^-function [29 . It has the following relation with the Weierstrassian a- 

function denoted by a w (u): <j w (u) oc e VlU a(u), rji = 7r 2 (i — ^Y^=i ) and q — e v ^ Tr . Consequently, 
our (^-function (|2.3[) is different from the Weierstrassian ^-function by an additional term ~2r\\u. 
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Here R 2 i(u) = Pi2Rf 2 (u)Pi 2 with P i2 being the usual permutation operator and U denotes 
the transposition in the i-th space. Here and below we adopt the standard notation: for 
any matrix A e End(C n ), Aj is an operator embedded in the tensor space C n <S> C n <S> • • •, 
which acts as A on the j-th space and as an identity on the other factor spaces; Rij{u) is an 
embedding operator of R-matrix in the tensor space, which acts as an identity on the factor 
spaces except for the i-th and j-th ones. The quasi-classical properties (12. 7p of the R-matrix 
enables one to introduce the associated classical Z n r-matrix r(u) [30] 

R B (u) = id + w r(u) + 0(w 2 ), when w — ► 0, 

1 — n> w \ cr/ (0)cx Q ,(u) 1 ... d , . . . 

An) = — c(«)+ e dR^ 7 -®^ • "<°> = < 2 - 8 > 



In the above equation, the elliptic C-function is defined in f |2.3j) . 

One introduces the "row-to-row" monodromy matrix T(u) [5], which is an n x n matrix 
with elements being operators acting on (C 1 )®^ 

T(u) = Rg(u + z 1 )R$ 2 (u + z 2 )---R$ N {u + z N ). (2.9) 

Here {zi\i — 1, • • • , iV} are arbitrary free complex parameters which are usually called inho- 
mogeneous parameters. With the help of the QYBE (II. ip . one can show that T(u) satisfies 
the so-called "RLL" relation 

Rg(u - v^u^v) = r 2 (v)Ti(«)i2f 2 (« - v). (2.10) 

An integrable open chain can be constructed as follows [6]. Let us introduce a pair of 
K- matrices K~(u) and K + (u). The former satisfies the RE 

Rf 2 { Ul - u 2 )Ki(u 1 )R% 1 (u 1 + u 2 )K 2 {u 2 ) 

= K 2 (u 2 )R? 2 { Ul + u 2 )Ki{ux)R% x {ui - u 2 ), (2.11) 

and the latter satisfies the dual RE 

Rf 2 {u 2 - u^K^u^R^i-Ui -u 2 - nw)K^{u 2 ) 

= K+(u 2 )R? 2 (- Ul -u 2 - nw)K+{ Ul )Rg{u 2 - u x ). (2.12) 

For the models with open boundaries, instead of the standard "row-to-row" monodromy 
matrix T(u) (12.91) . one needs the "double- row" monodromy matrix T(u) 

T(u) = T(u)K-(u)T~ l (-u). (2.13) 



Using (12.101) and (12. lip , one can prove that T(u) satisfies 

Rg( Ul - n 2 )T 1 ( Wl ) J Rf 1 ( Wl + u 2 )T 2 (u 2 ) = T 2 (u 2 )Rf 2 ( Ul + u^T^R^ - u 2 ). (2.14) 

Then the double-row transfer matrix of the inhomogeneous Z n Belavin model with open 
boundary is given by 

t(u) = tr(K + (u)T(u)). (2.15) 

The commutativity of the transfer matrices 

[r(«),r(v)] = 0, (2.16) 

follows as a consequence of (Il.ip - (I2.6I) and (I2.lip - fl2.12p . This ensures the integrability of 
the inhomogeneous Z n Belavin model with open boundary. 

In this paper, we consider a non-diagonal K-matrix K~(u) which is a solution to the RE 
(12. lip associated with the Z n Belavin R-matrix [31] 



The corresponding dual K-matrix K + (u) which is a solution to the dual RE (I2.12p has been 
obtained in [32J. With a particular choice of the free boundary parameters with respect to 
K~(u), we introduce the corresponding dual K-matrix K + (u) 

x0i:u,(-«)^UtM- (2.i8) 

In (I2.17P and ( 12. 18ft . 0, 0, are intertwiners which will be specified in section 4. We consider 
the generic {A,} such that \ ^ Xj (modulo Z + rZ) for i ^ j. This condition is necessary for 
the non-singularity of K ± (u). It is convenient to introduce a vector A = ^ i€i associated 
with the boundary parameters {Aj}, where {e,, i — 1, • • • ,ra} is the orthonormal basis of the 
vector space C n such that ( = 5ij. Moreover, in the following we always assume that 
the suffix index of the parameter Aj takes value in Z n cyclic group, namely, 

\±n = A i; i = l,...,n. (2.19) 



6 



3 ^n-i SOS R-matrix and face-vertex correspondence 

The Ai-i simple roots {c^} can be expressed in terms of the orthonormal basis {e^} as: 

ai = €i- e i+ i, i = 1, • • -,n - 1, 
and the fundamental weights {Aj | z = 1, • • • , n — 1} satisfying (Aj, otj) = 5ij are given by 

i . n 

k=l k=l 

Set 

I " 

/ e,; — e, e=— ^^Cfc, i = l,---,n, then ^^2 = 0. (3.1' 



n 

fc=i i=i 



For each dominant weight A = Y^=i a «Aj, a % ^ ^ + (the set of non- negative integers), there 
exists an irreducible highest weight finite-dimensional representation Va of A n _x with the 
highest vector |A). For example the fundamental vector representation is V\ ± . 

Let f) be the Cartan subalgebra of A n _i and f)* be its dual. A finite dimensional diagonal- 
izable f)-module is a complex finite dimensional vector space W with a weight decomposition 
W = ®fi^*W\pi], so that f) acts on W\p] by xv = /i(x)v, for any x G f), f G 
For example, the fundamental vector representation V\ 1 = C n , the non- zero weight spaces 
W[i] = Cei, z = 1, 

For a generic m G C n , define 

m i = (m,e i ), = - = (m, e* - e,) , z,j = l,---,n. (3.2) 

Let i?(u,m) G £nd(C n <g> C n ) be the R-matrix of the A^L X SOS model [53], 



i2(u, m) = ^ R%{u, m)E u ®E U + J2 m ) E u ® % + ^j( u ' ® #v } - ( 3 - 3 ) 

i=l i^j 

where E^ is the matrix with elements {Eij) l k = 5jk$u- The coefficient functions are 

4< (M , m) = ^M^ , (3.5) 
7 a(tt + w)a(mij) 
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and rriij are defined in (13. 2p . The R-matrix satisfies the dynamical (modified) QYBE 

Ri2( u i — u 2 ,m — wh^)Ri 3 (ui — U3, m)R 2 3(u2 — u 3l m — wh^ 1 ') 

= ^23(^2 - «3, m)R 13 (u l —u s ,m — whP^R^iui - u 2 , m), (3.6) 

and the quasi-classical property 

R(u,m)\ w -+o = id. (3.7) 

We adopt the notation: R^u, m — wm 3 ') acts on a tensor ®t> 2 <8> v% as i?(tt, m — wfi)®id if 
t>3 G W^f/i]. Moreover, the R-matrix satisfies the unitarity and the modified crossing-unitarity 
relation. 

Let us introduce n intertwiner vectors which are n-component column vectors <p m ,m-wj{u) 
labelled by j (j — 1, . . . , n). The k-th element of 4>m, m - W j{u) is given by 

e {k \u + nmj). (3.8) 

We remark that the n intertwiner vectors Ym,m-wj{u) are linearly independent for a generic 
m G C n . 

Using the intertwining vector, one derives the following face-vertex correspondence rela- 
tion [33] 

-^12(^1 ~ M a) y m,m—wi («l)® Ym— wi,m— w(i+j) 

= #K - U2, m)% r m -wi, m -w(i+k)( u i) ® </V„K^). (3.9) 

Then the QYBE (II. ip of the Z n Belavin's R-matrix R B {u) is equivalent to the dynamical 
Yang-Baxter equation (I3.6P of the A^ 1 ^ SOS R-matrix R(u,m). For a generic m, we may 
introduce other types of intertwiners 0, satisfying the conditions, 



6 (fc) 

m,m—wj\ / 



k=l 

fe=l 

from which one can derive the relations 



Cjm-wAM = <W» ( 3 - 10 ) 

fc=l 



TV -(u)0 (j) = (3.12) 



/x=l 
n 



^m+iufcmO^) ^m+w)/i,m( n ) — (3.13) 



With the help of (I3.10I) - (I3.13I) . we obtain the following relations from the face- vertex corre- 
spondence relation (13.91) : 

{<t> m +wk,m( U l) ® id ) #fa(«l - U 2) ( id ® ^m+uy.mK)) 

= XI ~~ M2 ' ® ^ m + W (fc+0, m+U) fe( U 2), (3.14) 



tm+w^mM ® ^m+^Cfe+O.m+^feK) J #i 2 («l ~ ^) 

= ^ jR(wi - M 2 , ^m+io(?+i),m+iuj(Wl) ® 4>m+wj,m {va) , (3.15) 



id ® ^m,m-«;f( M 2)J #12 («1 ~ "2) ( 0m,TO-m(wi) ® id 

= X #K - M 2, m)g TO _ w f )m _ TO (ifc + r)(wi) ® (f>m-wi,m-w{i+j)(u 2 ), (3.16) 

= ^ R(U! - U 2 , <Pm,m-wi{Ul) ® 4>m- w i, m -w(i+j)(u 2 ). (3.17) 

The face- vertex correspondence relations (I3.9P and (I3.14p -( l3.17p will play an important role 
in translating formulas in the "vertex form" into those in the "face form" . 

Corresponding to the vertex type K- matrices (I2.17P and (12.181) . one introduces the fol- 
lowing face type K-matrices /C and JC 



K(A|u){ x-MK(uy t( f>U x _ wi (-u), (3.18) 

s,t 

mri =E^S?a-^(-«) ^(m,, (3-i9) 



s,t 



Through straightforward calculations, we find the face type K-matrices simultaneously have 
diagonal forme 

K{X\u){ = 51 fc(A|«; Oi, £(A|i# = 8f k(\\u) u (3.20) 
where functions k(\\u; fc(A|w)i are given by 

fc(A|M;Ci = -7T v (3.21) 

a(Ai + t, +u) 

fc(A «)i = <^ [ [ 7T > > (3.22) 



3 As will be seen below (see (|4.27l) ). the spectral parameter u and the boundary parameter £ of the reduced 
double-row monodromy matrices constructed from IC(X\u) will be shifted in each step of the nested Bcthe 
Ansatz procedure |21j . Therefore, it is convenient to specify the dependence on the boundary parameter £ 
of JC(\\u) in addition to the spectral parameter u. 
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Moreover, one can check that the matrices K,(X\u) and K,(X\u) satisfy the SOS type reflection 
equation and its dual, respectively [32J. Although the K-matrices K ± (u) given by (I2.17P and 
(I2.18P are generally non-diagonal (in the vertex form), after the face- vertex transformations 
(I3.18P and (I3.19p . the face type counterparts fC(X\u) and JC(X\u) become simultaneously 
diagonal. This fact enables the authors to apply the generalized algebraic Bethe Ansatz 
method developed in [21] for SOS type integrable models to diagonalize the transfer matrix 
t{u) ([215]). 

4 Algebraic Bethe Ansatz 

By means of (13.121) . (13.131) . (I3.19P and (I3.20p . the transfer matrix r(u) (12.151) can be recast 
into the following face type form: 

t{u) = tr(K + (u)T(u)) 

= '^tr{K + {u)(j)x- w ^-u),\- w f l {u)^x-w(^ 

fj,,!/ 

= ^ (j)x,\- W fjj{-u)K + (u) (j)x- w (p,-u),X-wp,(u) (f)X-w(p-u),X-wjl(u) T{u)(j)x,X-wfi{- U ) 

= j2^Mur u nx\u^ = j2HMu),nx\u^. (4.i) 

Here we have introduced the face-type double-row monodromy matrix T(X\u), 

T(X\u)^ = <f>X-w(fi-i>),X-wfi(u) T(u) (px,x-wft(-u) 

= E*V,),a_ a HT(^<U A (-^)- ( 4 -2) 

This face-type double-row monodromy matrix can be expressed in terms of the face type 
R-matrix R(u, X) ([S3]) and the K-matrix K(X\u) (ETTgjl [2]]. Moreover from fl2TT4D . (1331) and 
(I3.13P one may derive the following exchange relations among T(m\u)^: 

x R( Ul + u 2 , m)fX r ( m + w (h + H)\U2)% 

*1,*2 jl ,31 

x T(m + w{h + 2a)k)£ - «2, m)gg. (4.3) 
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Following [21] and motivated by our recent work [20] for the open XXZ chain, we introduce 
n sets of operators {A^ s \ B^ S \C^ S ' labelled by s — 1, . . . , n, as follows: 

^(m|?i)=T(m|^, igf ) (m|^) = T(m| " )l , ^(ml«) = T(w| f , i ^ s, (4.4) 

(J [7 f lis) (T [Jflgi ) 

V^\(m\u) = ^ mjs - 6 * jW "> {T(m\u)i - 5lR(2u,m + ws) j s ; A {s) (m\u)} , i,j ^ s. (4.5) 

Some remarks are in order. Among the n sets of operators, the first set (corresponding to 
s — 1) is the very one which was used in [2T1 [23] to construct the algebraic Bethe Ansatz. 
Such algebraic Bethe Ansatz based on the first set of operators only gives rise to the first 
set of eigenvalues and associated Bethe Ansatz equations. In order to find the complete sets 
of eigenvalues, we find that the whole n sets of operators {A^ s \ B^ S \C^ S \V^\ s — 1, . . . , n} 
are needed. 

After tedious calculations analogous to those in [21], we have found the commutation 
relations among A^ s \m\u)^ T>^ s \m\u) and B^(m\u) from (14. 3p . Here we give those which 
are relevant for our purpose 

A {s \m\u)B\ s \m+w(i - s)\v) 

cr(u + v)a(u — v — w) „( s \ , .„ . ./„■>, .„ . 

= —. ' K . . {B\ >(m+w(i - s)\v)A {s >{m+w{i - s)\u) 

a(u + v + w)a{u — v) 

a(w)a(2v) a(u—v—m si +w)„( s ), /A ASI . ./_%, AS , s 

- a[ u- v l { 2rL) a( m „ - w) ' W^-s)M^>(m+^-*)H 

a(w) y a{u+v+m< ,, + 2w) B M {m+w{& _ i)]u)v (.r {m+w{i _ 3){v) 

< J rr m -i- id\ 



a(u+v+w) a[m as + w) 



i ± s, (4.6) 



V^imlujB^ (m + w{] - s)\v) 

a{u — v + w)a{u + v + 2w) 



x 



a{u — v)a{u + v + w) 

R(u + v + w,m - wi) K a2 P ^R(u — v , m + wj)^ t 

x Bf 2 (m + w{k + & - i - s)\v)V {s) 2 (m + w(j-s)\u)\ 

a(w)a(2u + 2w) j ^ a(u - v + m sa - w) kp 

> t r R(2u + w , m — wi) „7 



cr(u — v)a(2u + w) I ^— ' cr(m sa — w 



x ^ s) (ra + w(A; + /3- 2- s)|w)£> (s) J(m + s)|v)J 
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a(w)a(2v)a(2u + 2u>) 



a(u + v + w)a{2v + w)a(2u + w) 
a(u + v + m S j^ 




R(2u + w,m — wi) 



a(m S j — w) 

x B$(m + w(k + a — i— s)|-u)^l (s) (m + w(j - s)\v) 
i,j,k i- s, 



(4.7) 



B\ s \m + w{i — s)\u)Bj S) \m + w(i + j — 2s)\v) 

= ^R{u-v,m-2wsf j iBf(m+w0-s)\v)B^\m+w{a+^ 

a,f3=£s 

hi ^ s. 



(*)/ 



For the special case of s = 1, the above commutation relations fl4.6p - fl4.8l) recover those in 



In order to apply the algebraic Bethe Ansatz method, in addition to the relevant com- 
mutation relations fl4.6l) - fl4.8p . one needs to construct a reference state associated with each 
s, which is the common eigenstate of the operators ^4^, T>^\ and is annihilated by the 
operators Q . In contrast to the trigonometric and rational cases with diagonal K ± {u) [B], 
the usual highest-weight state 



1 1 ^ 







® ■ • • ® 


V : J 





is no longer the pseudo-vacuum state. However, after the face- vertex transformations (13.181) 
and (13.191) . the face type K-matrices JC(X\u) and K,(X\u) simultaneously become diagonal. 
This suggests that one can translate the Z n Belavin model with non-diagonal K-matrices 
into the corresponding SOS model with diagonal K-matrices K,(X\u) and K,(\\u) given by 
(I3.18p - ()3.19p . Then one can construct the pseudo- vacuum in the "face language" and use 
the algebraic Bethe Ansatz method to diagonalize the transfer matrix. 

One of the reference states corresponding to the case of s = 1 (or the first one of the n 
reference states ( 14. 9 p below) was found in [21] and yields the first set of eigenvalues of the 
transfer matrix. Here, we give the complete n reference states {|f2^(A))|s = 1, . . . , n}. For 
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each s, we propose 

|fi (s) (A)) = 0A-(iV-l)«)8,A-JVti>s(-«l) ® 0A-(JV-2)«,s,A-(7V-l)«,s(-2 ; 2) ' ' ' ® 0A,A-™s( _;j iv) ■ (4.9) 

These states (s = 1, . . . , n) depend on the boundary parameters {Aj} and the inhomogeneous 
parameters {zj}, but not on the boundary parameters £ and £. We find that the states given 
by (14.91) are exactly the reference states in the following sense, 



V^){\-Nws\u)\rt s \\)) = 3/«(„)fc(A|u+^;£-^ 



x ■ 



^ S )(A- A%s|w)|ft (s) (A)> = £;(AKOs|fi (s) (A)), (4.10) 

ft , ^ + ^~^ A |nW(A)>, 

ijVs (4.11) 

C, (s) (A- iVws|u)|fi w (A)) = 0, z^s (4.12) 

4 s) (A-^s|u)|fi (s) (A)) ^ 0, i^s. (4.13) 



Here f^ s \u) is given by 



r (•)/-..■, _ °"( 2 «)o-(A s + u + w + f) 



/S(M) <t(2u + u;)<7(A 8 + « + (4 ' 14) 



In order to apply the algebraic Bethe Ansatz method to diagonalize the transfer matrix, 
we need to assume N = n x I with / being a positive integer [2T]. For convenience, let us 
introduce a set of integers: 

AT. = (n-i)xl, i = 0,l,---,ra-l, (4.15) 

and n ( n ~ 1 ) / complex parameters {i>^| = 1, 2, • • • , A^ +1 , i = 0, 1, • • • , n — 2}. As in the usual 
nested Bethe Ansatz method, the parameters {i>^} will be used to specify the eigenvectors 
of the corresponding reduced transfer matrices. They will be constrained later by Bethe 
Ansatz equations. For convenience, we adopt the following convention: 

v k = vf\ fc = l,2,...,iV 1 . (4.16) 



4 Such states played an important role in constructing extra centers of the elliptic algebra at roots of unity 

m. 
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We will seek the common eigenvectors (i.e. the so-called Bethe states) of the transfer matrix 
by acting the creation operators on the reference state |f2^(A)) 

\v lr --,v Nl )M = ^ 1 ' l2 '---' i ^Si 1 s) (A+n;^-s)|^ 1 )Si 2 s) (A+u;(z 1 +2 2 -2s)|^ 2 ) 
h,---,iN 1 ¥ :s 

Ni 

x • • • B« (A + w h ~ wN x s\v Nl )\rt s \\)). (4.17) 
fc=i 

The indices in the above equation should satisfy the following condition: the number of 
i k = j, denoted by #(j), is I, i.e. 

#0') = f. 3^s. (4.18) 
Then (13. ip and the above restriction (14.181) imply 

Ni n 

\ + w^2i k - wNJ = X + wl^j-wNxs = \-w(l + N^s = A - wNs, (4.19) 

k=l j^s 

which is crucial for the diagonalization of the transfer matrix in the remaining part of the 
paper. 

With the help of (14. ip . (I4.4p and (14. 5 p we rewrite the transfer matrix (12.151) in terms of 
the operators A^ and X>W 

n 

r(«) = ^k{\\u) v T{\\u)l 

= k{X\u) s A {s) {X\u) + ~ k ( x \ u )i T W u ) i i 
= ~k(X\u) s A (s) {X\u) + J2HMu)iR{2u,X + ws) i s s i A < - s \X\u) 
+ ^fc(A|u) i (t(X\u)\ - R(2u,X + wsy°A is \X\u) 

n 

= Y.HM^iR^X + wsY^A^iXlu) 
1=1 

+ ^fcW(Aiu+|) ^^ 

V^8 



= «?>(«) A s \X\u) + ^(Ak + ^)*V^\{X\u). (4.20) 
Here we have used (14.51) and introduced the function otpiu), 



J2HMu)iR( 2 u,X + ws)% (4.21) 



i=i 
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and the reduced K-matrix K^(X\u) with elements given by 

icU(\\u)i = Sf~k<»(X\uk, i,j^s (4.22) 

mMu)i = { n %=fU t'!r M+ S!!| - 1 * a (4 - 23) 

Using the technique developed in [2T], after tedious calculations, we find that with the 
coefficients F ll ' Z2 '"'' iN i in (14.171) properly chosen, the Bethe state |t>i, • • • , f _/Vi ) becomes the 
eigenstate of the transfer matrix (12.151) . 



T 



{u)\ Vl , ■ ■ • , ^) (S) = A s (u; e, K»K ■ • ■ , ^) (s) , (4.24) 



provided that the parameters {v® \ k = 1, 2, • • • , iVi+i, z = 0, 1, • • • , n— 2} satisfy the following 
Bethe Ansatz equations: 

(1) a(2vj + w) j+ a(vj + v k )a(vj - v k - w) 



cr(2uj + 2w) , A f-_ 1 cr(fj + u k + 2w)a(vj - v k + w) 



Mi 

N 



= n gEl±£^^L A M fe + «. { _ » {„«}), ( 4 . 25) 
Hj- cr(vj + z k + w)a(vj - z k + w) 2 2 

3 ' <r{7vf + 2w\} 3 l =l a{vf + + 2w)a{vf - vf + u,) 

k=1 a{v)'+zl'+w)a{v)'-zl'+w) z * 

i = l,---,n-2. (4.26) 

Here {/3i^(tt)|« = l,...,n — 1} are functions given in Appendix A, {A s (u;£,{v^})\i = 
1, . . . , n — 1} are the eigenvalues (given in Appendix A) of the reduced transfer matrices in 
the nested Bethe Ansatz process, and the reduced boundary parameters and inhomo- 

geneous parameters {z k } are given by 

^) = £0-|, + < = ,...,n-2. (4.27) 

In the above we have adopted the convention: £ = zj® = z k . The corresponding 
eigenvalue A s (u; £, {v k }) is given by 

c(A s + § + u) jj^JL cr(w + w fc + w)<t(m - ujfe) 
15 



^a(2u)a(X s + u + w + £) -i-j- <t(w — v k + w)a{u + t> fc + 2w) 
cr(2u + w)<t(\ s + u + £) J-jL <t(m — v k )a(u + Vk + w) 

xTT ? g( " + ^f- zt) , Ai 1) (^+-;e--,K (1) })- (4.28) 

It is easy to check that the first set of eigenvalues Ai(u; ^, {v^}) and the corresponding 
Bethe Ansatz equations are exactly those found in [21\. However, the rest n — 1 sets of 
eigenvalues {A s (u; £, {wfc})|s = 2, . . . , n — 1} and the associated Bethe Ansatz equations are 
new ones. For the special case of n = 2, which corresponds to the open XYZ spin chain, 
we find that the two sets of eigenvalues {A s (u; £, = 1, 2}, after rescaling of an overall 

factor due to different normalizations of the R- and K-matrices, recover those in [35] obtained 
by directly solving the T-Q relation (or functional Bethe Ansatz). As shown in [35], these two 
sets of eigenvalues give the complete spectrum of the transfer matrix of the open XYZ spin 
chain. As a consequence, the corresponding two sets of Bethe states {\vi, ■ ■ ■ ,v Nl )^\s = 1,2} 
together constitute the complete eigenstates of the transfer matrix of the open XYZ model. 
Therefore, it is expected that the n sets of eigenvalues {A s (u; £, {vfc})|s = 1, . . . , n} (14.281) 
[resp. Bethe states {\vi, • • • , v^p"'^ = 1, . . . , n} (I4.17P and (I4.25p - (l4.26p ] together give rise 
to the complete spectrum [resp. the complete eigenstates] of the transfer matrix (I2.15P of 
the open Z n Belavin model. 

5 Result for the associated Gaudin model 

As will be seen from the definitions of the intertwiners (13. 8p . (13.101) and (13. lip , specialized to 
m = A, (j)\ t \- w i{u) and 4>\,\- w i(u) do not depend on w while 0a,a-iot(w) does. Consequently, 
the K- matrix K~(u) does not depend on the crossing parameter w, but K + (u) does. So we 
use the convention: 

K(u) = lim K~(u) = K~(u). (5.1) 

w— >0 

We further assume that the parameter £ has the following behaviojfl as w —>■ 0, 

£ = £ + w S + 0(w 2 ), (5.2) 

5 In [TSJ [23l [20] , a special case of £ = £ was studied. The generalization to the case with nonvanishing S 
is straightforward. 
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with a constant 5. It implies that 



lim £ = £ . 

u)->0 

Then the K-matrices satisfy the following relation 

lim {K + (u) K~(u)} = \im{K + (u)}K(u) = id. (5.3) 

w^O m— »0 

Let us introduce the elliptic Gaudin operators {Hj\j = 1, 2, • • • , A^} associated with the 
inhomogeneous Z n Belavin model with open boundaries specified by the generic K-matrices 
OTTj) and fT2~T8|) : 

N ( N } 

Hj = Tj(zj) + ^2 r kj( z i ~ z k) + K j l {zj) \ ^2 r i k ^ z i + Zk ) \ K A z i)' ( 5 - 4 ) 

k^j { k+j ) 

where Yj{u) = £{Kj(u)}\ w=0 Kj(u), j = l,--.,N, with Kj (it) = tr {K+(u)R^(2u)P oj } . 
Here {zj} are the inhomogeneous parameters of the inhomogeneous Z n Belavin model and 
r(u) is given by (12.81) . For a generic choice of the boundary parameters {Ai, • ■ ■ ,X n , £}, 
Tj(u) is a non-diagonal matrix. 

Following [521 EZ], the elliptic Gaudin operators (15. 4p are obtained by expanding the 
double-row transfer matrix ( 12. 15ft at the point u = Zj around w — 0: 

t(zj) = T(zj)\ w=0 + wHj + O(w 2 ), j = l,---,N, (5.5) 
d 

Hj = -^T(zj)\ w=Q . (5.6) 

The relations (12.71) and (I5.3p imply that the first term r(zj)\ w =o in the expansion (15. 5p is 
equal to an identity, namely, 

T(zj)\ w=0 = id. (5.7) 

Then the commutativity of the transfer matrices {t(zj)} (12.161) for a generic w implies 

[H J ,H k }=0, i,j = l,.--,N. (5.8) 

Thus the elliptic Gaudin system defined by ( 15. 4p is integrable. Moreover, the relation (15.61) 
between {Hj} and {t(zj)} and the fact that the first term on the r.h.s. of (15. 5p is iden- 
tity operator enable us to extract the eigenstates of the elliptic Gaudin operators and the 
corresponding eigenvalues from the results obtained in the previous section. 
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Using (lOTj) . (pOp , (TOD and (CO]) -([Oh the Bethe Ansatz equations P5| and (14^261) 
become, respectively, 

^ ( ^) } ^(2f +tt>) || + vf)a(yf - vf - w) 



a(2vf + 2w) k ^ti a{vf + 4 ] + 2w)a(vf - vf + w) 



' ' / (*) i i 3w\ ( (») (»-l) i 



k=l u \ v j ^ v k ^ 2> u y u j u k ^ 2> 

x jf +4 i+1) +iH-f-e i) -f: 



tv(,f +4 i+1) + fV(^-4 m) + f)' 

z = 0,---,n-3, (5.9) 

/ (n-2) . \ iV„_i / (n-2) . (n-2)x , (n-2) (n-2) x 

cr(2t>j ; +u>) -p-j- £r(uj ; + ^ )er(vj - v\ ' — w) 



mr*) r u :: ii 



/r, (n-2) , r, \ 11 / (n-2) . (n-2) . x / (n-2) (n-2) . x 

a(2v) >+2w) k , jk= 1 a{v) > +v k ' +2w)a{v) >-v k +w) 



" ^(a s _! + e - vf- 2) - W )a{\ s ^ + ^-u + vf~ 2) + f ) 

x l\ M n - 2) + -t 3) + »M«r> - 4 n " 3) + f ) ' 1 ' j 

Here we have used the convention: vf 1 ^ = Zk, k = 1, • • • , N. The quasi-classical property 
( 13. 7p of R(u, m), (1A.1I) and (|A.2[) lead to the following relations 

/f +1 )(u,£,0) = l, |-/3i m )(n,e,0) = 0, i = 0,-..,n-2. (5.11) 

Noticing the restriction (15.21) . one may introduce some functions {7s (w)} associated with 

Using ( 15.51) . we can extract n sets of eigenvalues {h^\s — 1, . . . , n) (resp. the corresponding 
Bethe Ansatz equations) of the Gaudin operators Hj ( 15.41) from the expansion around w = 
for the first order of w of the eigenvalues (14.281) of the transfer matrix r(u = Zj) (resp. the 
Bethe Ansatz equations (I5.9P and (15.101) ). Finally, the eigenvalues of the Z n elliptic Gaudin 
operators are 

hf = iPizj) + C(A S + e + Zj ) - {C(% + Xk) + C(% - Xk)} , (5.13) 

fe=i 
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= — I w)\ w=0 + 6 — F P { : +1 \u, 1 0)1^, i = 0, • • • , n - 2. (5.12) 



where ^-function is defined in (12.31) . The ra( -" 2 — parameters {x k ^\ k = 1,2,- ■■ ,N i+ i, i = 
0, 1, • • • , n — 2} (including X 3jS X — X k , k = 1, - ■ ■ ,Ni) are determined by the following 
Bethe Ansatz equations 

N i+1 
Ni 

= 7r 2) (4 i} ) - E {c(*f + 4 i_1) ) + ccf - 
fe=i 

ATi+2 

+c(Wi+e+xf)-E {c(4 J) +4 m) )+c(4 J) -4 m) )} , 

fc=i 

i = 0, 3, (5.14) 

7 (n-l) (x (n- 2 ) ) _ c(2 ^ 2)) _ 2 £ {c(xf- a )+X^) + C(x^-xi"- 2 ))} 

= c(A^ 1+ e+xf- 2) )+(^ - *) c(A s _ 1+ e-xf- 2) ) 

- e {c(4 n - 2) +4 n - 3) )+c(xr 2) -xr 3 ))} . (5.i5) 
fe=i 

Here we have used the convention: x k — z^, k — 1, • • • , N in (I5.14p . Then the n sets 
of eigenvalues {/ij |s = l,...,n} given by (15. 13[) - (15. 15[) (c.f. [23]) together constitute the 
complete spectrum of the Gaudin operators Hj (15.41) . 

6 Conclusions 

We have discovered the multiple reference state structure of the Z n Belavin model with 
boundaries specified by the non-diagonal K- matrices K ± (u), (12.171) and (1 2 . 1 8 D . It is found 
that there exist n reference states {|f2( s )(A)) \s = 1, . . . , n} (14.91) . which lead to n sets of Bethe 
states \vi, ■ • - , v^)^ (14.171) . These Bethe states give rise to n sets of Bethe Ansatz equations 
(I4.25p - (l4.26p and eigenvalues (14.281) . labelled by s = 1, . . .,n. The fist set of them, which 
corresponds to the s = 1 case, gives the results found in [21]. It is expected that these n sets 
of eigenvalues {A s (u; £, {vk}) \s = 1, . . . , n} (14.281) together give rise to the complete spectrum 
of the transfer matrix r(n) (12.151) for the Z n Belavin model with generic boundaries. In the 
quasi-classical limit (i.e. w — > 0), the resulting n sets of eigenvalues {h^\s = 1, . . . ,n} given 
by (I5.13I) - (I5.15I) together constitute the complete spectrum of the Gaudin operators Hj (15 .4p . 
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Taking the scaling limit [38] of our general results, for the special n = 2 case, we recover 
the results obtained in [20J for the open XXZ spin chain. It is believed that such structure 
of multiple reference states also exists for the open A^li trigonometric vertex model studied 
in 
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Appendix A: Definitions of the a-, (3-, A- functions 

In this appendix, we give the definitions of the functions a^\u), $\u) and A«( u; £,{««}), 
which appeared in the expressions of the eigenvalues and the Bethe Ansatz equations (14.251) - 
(14281) . 

In order to carry out the nested Bethe Ansatz for the Z n Belavin model with the generic 
open boundary conditions, one needs to introduce a set of reduced K-matrices K,^\\\u) 
which include the original one K,(\\u) = JC^(\\u): 

£ {r \\\u){ = 5i¥ r \X\u)i, i,j = r + l,---,n, r = 0,---,n-l, 

Moreover we introduce a set of functions {a^ r \u)\r = 1, ■ ■ • , n — 1} related to the reduced 
K-matrices KM\\\u) 



^R{2u 1 \ + wr)%U r - l \\\u) i , r = l,---,n- 1, (A.l) 



and an associated set of functions {/3^'(u, £,w)\i = 1, . . . ,n — 1} 

^(u,^)=P^\u) = a^\u) " (A,;+ ;^"!~ fW - v < = 0,-.,n-2. (A.2) 

cr(A i+ i + 4 + u + w - pv) 

In the process of carrying out the nested Bethe Ansatz [21], one needs to introduce a set 
of functions {AW(iz; {v^ — 0, . . . , n — 1} which correspond to the eigenvalues of the 
reduced transfer matrices. The functions {A^(u; £, {i^ })} are given by the following recur- 
rence relations 

a® {*$>}) = («) a{x r: + e( !n u + , w) n g(u + t Mm = - = $ 

fe=i 17 ( M + ^fc + w)a{u — Vf. ) 
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a{2u+w)a(\ i+1 +uH {i) ) Li <t(u-v®)<t(u+v®+w) 



k=l ^V" ^A; ) u \ u '^ u k 

xfl /'"^"w""!'' A" +1, ("+f^" ) -f.H' +1 '}). 

z = l,---,n-2, (A.3) 

o"(A n + 4 - « - f M\» + 6 } + tt) 

The reduced boundary parameters and inhomogeneous parameters {z^} are given by 

(l427j) . It is remarked that all the functions {a w (w)}, {/3 (i) (w)} and {A^(u; £, {v^})} are 
indeed the functions of the boundary parameters {Aj}. 

Since that the suffix index of the boundary parameters {Aj} takes value in Z n , one can 
introduce a Z n cyclic operator P (i.e. V n = id), which acts on the space of functions of 
{Ai, . . . , A n }. On any function f(Xi, . . . , A n ) the action of the operator V is given by 

T 3 (/(Ai, • • • , A n _i 

j A n )) — /(Ai+i, . . . , A n , A n+ i) — /(A2, . . . , A n , Ai). (A. 5) 

Then we introduce the following a-, /3-,A-functions: 

a W( u ) = p 5 - 1 (a (i) (u)), s = l,...n, z = l,...n-l, (A.6) 
/?«(«) = P^ 1 (/?«(«)), s = l,...n, i = l,...n-l, (A.7) 
WteiMPY) = V s - 1 ^{u^,^}), s = l,...n, z = l,...n-l. (A.8) 

Direct calculation shows that the two definitions of a« , (14.211) and (1A.6I) . coincide with each 
other. 
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